Surface effects in multiband superconductors. Application to MgB2. 
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Metals with many bands at the Fermi level can have different band dependent gaps in the su- 
perconducting state. The absence of translational symmetry at an interface can induce interband 
scattering and modify the superconducting properties. We dicuss the relevance of these effects to 
recent experiments in MgB2. 
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I. INTRODUCTION 

Recent experiments [Q report the existence of super- 
conductivity at nearly 40K in MgB2. Its origin is not 
completely elucidated. The material shows a pronounced 
isotope effect Q , and the density of states is well approx- 
imated by the BCS theory Tunneling experiments 
suggest that the superconducting properties at the sur- 
face of the material differ from the expected bulk behav- 
ior [^]. In addition, photoemission results |5| suggest 
the existence of an s-like gap. A, such that A < iksTc. 
A possible explanation of this result is that the measured 
A is the average of different gaps. The present work is 
motivated by the persistent discrepancy between the gap 
values measured in different experiments, and, particu- 
larly, the excellent fit to a BCS gap too low to explain 
the value of the critical temperature observed in tunnel- 
ing experiments reported in Q . 

Band structure calculations suggest that there are, at 
least, two types of bands at the Fermi surface: a hole 
band, built up of boron a orbitals, with a weak disper- 
sion in the direction perpendicular to the boron planes, 
and a broader band, built up mainly of tt boron orbitals, 
which shows a significant dispersion in the direction per- 
pendicular to the boron planes. Theoretical arguments 
favor, as the origin of the superconductivity, the hole 
like a band or the tt band |jl^. The existence of 

two bands with different physical properties is assumed in 
other models for the superconducting properties of MgB2 
It has been argued that the upper critical field 
can be best modelled if the superconducting properties 
depend on the specific band at the Fermi level |jl^. On 
general grounds, it is reasonable to assume that the a and 
TT bands in MgB2 will have different contributions to the 
superconducting phases, and that the superconducting 
gap needs not be the same in the two bands. 

The existence of many bands at the Fermi level, with 
very different physical properties, is probably a generic 
feature of intermetallic superconductors ||lj,|l^ . In these 
materials, it can be expected that the pairing interaction 
which gives rise to the superconductivity will depend on 
the details of each band. If this is the case, there is not 
a uniform gap at the Fermi level. The superconducting 
state resembles, in this respect, that of an anisotropic 
superconductor. The effects of interband scattering on 
the bulk properties of a superconductor with two different 



bands at the Fermi level was studied in |^ . 

Interband scattering can be induced by any defect 
which breaks the translational symmetry of the lattice, 
including lattice imperfections, phonons, and surfaces or 
grain boundaries. Thus, in materials with many bands 
crossing the Fermi level, one can expect the supercon- 
ducting properties near a surface to differ from those at 
the bulk. This effect should be more pronounced in sys- 
tems with many irregular interfaces, such as ceramic and 
granular materials. 

In the present work, we study pair breaking effects at 
the surfaces of many band superconductors. The main 
difference with the theory developed in |l^] is that, when 
interband scattering is restricted to the surface region, 
the superconducting gaps are inhomogeneous, because 
the value of the gaps within the bulk of the system is not 
affected by the presence of a surface. In order to take this 
effect into account, we develop a method which differs 
from the standard treatment of pair breaking effects in 
conventional [l7| or unconventional [l^ superconductors. 
We use the simplest possible model where non trivial ef- 
fects are expected, presented in the next section. Section 
III presents the main results, and the main conclusions 
which can be drawn from them are discussed in section 
IV. 



II. THE MODEL 



A. The hamiltonian. 



Superconductors with different bands at the Fermi 
level, where each of these bands have different pairing 
interactions have already been discussed in the literature 
[|l9|,0. The simplest model contains two bands, with 
two different densities of states and pairing interactions. 
We will consider the hamiltonian: 



H = Ho + Hint + Hph 
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Wpfc = E / rf'f^/(2)^L(^)^2,.(f) (1) 

Electrons within each band experience a different pair- 
ing interaction, 5^, leading to two superconducting gaps. 
The two bands are coupled by the interaction g' . Other- 
wise, in the absence of interband scattering the two gaps 
open at different temperatures, T^^i oc exp (— Wi/f/i), 
where Wi is the bandwidth, and is a cutoff related 
to the pairing mechanism. Specific heat measurements 
p]| seem to exclude this possibility in MgB2. We ne- 
glect intraband scattering, which docs not give rise to 
pair breaking effects, at least to lowest order j2^. Fi- 
nally, we assume a constant interband scattering term, 
localized near the surface. The function ]{z) is peaked 
near the surface, located at z = 0. The width of /(z) is of 
the order of a lattice spacing. With these restrictions, the 
model described in eq.(|^) includes six parameters, with 
dimensions of energy: the density of states at the Fermi 
level of each band, Niitp), the pairing interactions, gi 
and 5', and the scattering potential, V . 



B. Pair breaking effects. 




FIG. 1. Sketch of the lattice solved in the text. 

The BCS equations are solved using the standard 
equivalence of the attractive Hubbard model in a bipar- 
tite lattice to the repulsive Hubbard model [Q. The 
semiinfinite model is solved using transfer matrix tech- 
niques |2^], which are described in the Appendix. The 
gaps are calculated selfconsistently in a layer of m sites, 
which are coupled to an homogeneous chain where the 
gaps take the bulk values, where the gaps are defined as: 

We present results obtained with m = 32. This ideal- 
ized one dimensional model can be viewed as an approxi- 
mation to the inhomogeneous layered structure expected 
near the surface. Note that the BCS equations exclude 
the possibility of one dimensional fluctuations, so that 
the solutions to be discussed below do not show unphys- 
ical one dimensional features. 



The lack of translational symmetry induced by the 
surface makes it convenient to solve directly the the 
Bogoliubov-de Gennes equations derived from (|^). We 
use a discrete tight binding model for this purpose. We 
assume that each band can be described by a single or- 
bital per site, and that there are local attractive interac- 
tions with induce the pairing. The model reduces, in the 
absence of interband scattering, to two coupled negative 
U Hubbard models. In this basis, interband scattering 
can be included by allowing for hopping from one orbital 
to the other at any given lattice site. Our discretized 
model in real space becomes: 

H = Ho + 'Hint + Hph 
Ho = E ^A,^.,s^h],s + h.c. 

IS 

i 

Up, = E ^^,.,.^2,.,. (2) 

We assume that the lattice is a semiinfinite chain, and 
that interband scattering is restricted to the outermost 
site, as schematically depicted in Fig.[|l|]. 



III. RESULTS 

We use as unit of energy one of the hopping parame- 
ters in eq.(||), ti — \. We assume that band 1 represents 
the total effect of the two a bands in MgB2 . The density 
of states in the tt band is roughly one third that of the 
two a bands and we set t2 = Mi. We assume that the 
superconductivity is mostly due to the a band, and we 
choose the parameters Ui , U2 and U' so that the smallest 
gap is less than one half the largest one, at zero tempera- 
ture 1^. A reasonable combination is Ui = 0.5ti, U2 = 
and U' = 0.2ii so that Ai = 0.0212ii and A2 = 0.0079ii. 
The critical temperature in these units is Tc — O.Ollti. 
Using the experimental value Tc « 40K, the value of the 
unit of energy, ti is equal to 0.3 eV, and the density of 
states of the a band is 0.5 states/eV-cell, which has the 
right order of magnitude Note that, with these pa- 
rameters, the model is well into the BCS weak coupling 
regime. Finally, we must determine the strength of the 
interband scattering at the surface, V. The mixing of 
the boron tt and a bands is strongly suppressed in the 
bulk. This effect is weakened if the relative angles be- 
tween nearest neighbor B ions is modified at the surface, 
so that the a bands are not built up of sp^ orbitals only. 
The interband hopping associated with these deforma- 
tions will be a fraction of the bulk hopping terms. We 
take in the following V = OAti « 0.12eV, and restrict 
this interband hopping to the outermost layer. The gaps 
near the surface, at zero temperature, are plotted in the 
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inset in Fig.g. The gaps are almost uniform, suggest- 
ing that the influence of the interband scattering on the 
value of the gaps is very small. 



Interband scattering has a stronger influence on the 
density of states near the surface. The results for the 
outermost site are plotted in Fig.|^. The couphng be- 
tween the two bands shows in the existence of peaks in 
the density of states at the two gap positions, while, in 
the bulk, each band displays a single peak at the value 
of the corresponding gap. 
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T = 36 K 
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FIG. 3. Density of states at the surface at T = 36K. 
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FIG. 2. Density of states at the surface at zero tempera- 
ture. The inset shows the superconducting order parameters 
in each band as function of position from the surface. 



In some experiments, like point contact spectroscopy, 
the probe can be an additional source of interband scat- 
tering, at the position where the measurement is being 
made. In the presence of strong interband scattering at 
the surface, the perturbation in the densities of states is 
more pronounced, as shown in Fig.[^, calculated using 
V — ti K, 0.3eV. A single smeared gap will be observed 
in an experiment of this type. 
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FIG. 4. Density of states at the surface at T = OK, in the 
presence of strong interband scattering at the surface. 



The overall features in the density of states remain 
the same at relatively high temperatures, as shown in 
fig.|^, where the results at T=36K are shown. The scale 
at which the smallest gap closes is determined by the 
largest gap. 



It is interesting to note that, if the value of the gaps in 
the bulk were of opposite sign, a midgap state, induced 
by Andreev reflections, should arise. The problem con- 
sidered here maps onto that of a dimerized chain. If the 
value of the dimerization changes sign, a gap state always 
arises, as extensively discussed in connexion to solitons 
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in polyacetilene pst - 



VI. APPENDIX 



IV. CONCLUSIONS 



We have analyzed the effects of interband scattering at 
the surface, in a superconductor with a many bands at 
the Fermi level, and different pairing strengths for each 
band. 

Surface scattering, for reasonable values of the param- 
eters, is inneffective in changing the gaps near the sur- 
face. This is consistent with the fact that, in a conven- 
tional weak coupling BCS superconductor, the coherence 
length is much greater than the surface layer where strong 
scattering is expected. 

The influence of the surface scattering induces signif- 
icant changes in the density of states near the surface. 
The density of states of each band shows peaks at the 
positions of the gaps associated with the other bands. 
There is, however, always a minimum gap, below which 
the density of states vanishes. This gap, A^m, has a 
weak temperature dependence until close to the bulk crit- 
ical temperature. 

We have used a set of parameters appropiate for MgB2. 
We assume that there is a wide, delocalized band, and 
a narrower and more localized band, which determines 
mostly the superconducting properties. The wide band 
(derived from the tt orbitals) has the smallest gap, and 
is weakly influenced by the narrow band. The narrow 
band shows stronger features at both the small and large 
gaps. Tunneling experiments are highly sensitive to the 
delocalization of the wavefunction. Hence, it is possible 
that they measure the density of states associated to the 
wider band (tt orbitals). If this is the case, the observa- 
tions will show a single gap, well approximated by the 
BCS expression. This gap should have a weak temper- 
ature dependence, until temperatures comparable to the 
bulk critical temperature 

The strength of the interband scattering at the surface 
studied here can depend on the experimental setup, and 
it may be enhanced in some experiments. If that is the 
case, a single smeared gap will be observed, of magnitude 
comparable to 1.7 Tc |p7| . 

Finally, it is interesting to note that, if the gaps in 
the two bands were of opposite signs, as expected from 
electron-electron pairing mechanisms, a surface state 
near the center of the gap should appear. 



The problem defined by the hamiltonian in eq. (^) with 
the geometry shown in Fig.||l[] reduces to the calculation 
of the density of states in 2m sites, with gap values 
and connected by hopping terms which can be either i i , i2 
or V . The attractive Hubbard model, Ui < 0, can be 
mapped onto the repulsive Hubbard model, in a bipartite 
lattice, by the transformation: 

^iT "it 

cl ^ (-i)^da 

ca ^ (-i)Vlx (4) 

The superconducting order parameter, ^'^ = (c||c]^|), is 
mapped onto a staggered magnetization in the transverse 
direction. At half filling, there is an additional symme- 
try which allows us to rotate this magnetization to the 
z-axis. Then, '^i = ((— l)*(c||Ci| — c||Ca)). In this repre- 
sentation, the hamiltonian does not mix the spins, and it 
can be decomposed into two boxes, one for each spin di- 
rection. The problem is reduced to the calculation of the 
density of states in a tight binding chain with variable 
hoppings, ti i±i, and energy levels, ej, which are related 
to the local value of the gaps. The gaps must be deter- 
mined selfconsistency from the values of the ^iS. 
The fractions: 



9n,n±l{^) — tn,n±l 



Gn,n' (w) 



G 



?T,±l,n 



are independent of n', and satisfy: 



ffn,n±l(t^) = 

and: 

Gn,n{l^) = 



n , 71. ± 1 



(5) 



(6) 



(7) 

Thus, the problem can be solved by iteration from 
the boundaries, provided that one knows the values of 
g±m,±m±i- These values can be easily be calculated, if 
one assumes that the values of the e^'s are constant be- 
yond position m [^6| . Finally, the selfconsistency require- 
ment for the values of e;, i = 1, 2m must be satisfied. 
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